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Solutions to a class of differential systems that generalize the Halphen system 
are determined in terms of automorphic functions whose groups are commensurable 
with the modular group. These functions all uniformize Riemann surfaces of genus 
zero and have q-series with integral coefficients. Rational maps relating these func- 
tions are derived, implying subgroup relations between their automorphism groups, 
\Q . as well as symmetrization maps relating the associated differential systems. 

o 
o 

O , 1. Introduction 

00, 

Differential equations satisfied by modular functions have been studied since the time of 
Jacobi [J] . Such equations arise naturally in connection with second order Fuchsian differential 
operators whose monodromy groups coincide with the automorphism group of the function. 
Recall that a meromorphic function /(r) defined on an open, connected domain T> in the 
Riemann sphere is said to be automorphic with respect to a group 0/ of linear fractional 
transformations 

X: r: _£±» sr(T) , (:»W C) , (i.i) 



if the domain T> is (J5/ invariant and 

f(T(r)) = f(r), VT6 0/. (1.2) 

In this paper, the term modular will be applied to functions with Fuchsian automorphism 
groups of the first kind that are commensurable with the modular group PSL{2, Z) (i.e., 
whose intersection with the latter is of finite index in both). 

Using the terminology of [Fo], a "simple" automorphic function is one without essential 
singularities at ordinary points, whose fundamental region has a finite number of sides and 
which has a definite (finite or infinite) limiting value at any parabolic point (cusp). The 
following standard theorem gives the connection between linear second order equations and 
simple automorphic functions (cf. [Fo], Sec. 44, Theorem 15) 
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Theorem 1.1: If f(r) is a nonconstant, simple automorphic function, then the (multi- 
valued) inverse function r = r(f) can be expressed as the quotient of two solutions of a 
second order linear equation 

^ + R(f)y = o, (1-3) 

where R is an algebraic function of f . If f has a single first order pole in the fundamental 
region, then R is a rational functions of f . 

In fact, the solutions of (1.3) may be expressed, at least locally, as 

,=&$m. (i.4) 

( T ) 2 

Conversely, given any second order linear equation 

d 2 y „dy , nn n _v 

dp + W + Qy = ( 5) 

with rational coefficients P(f), Q(f), singular points (ai, . . . a n , oo), a basis of solutions (yi, 1/2) 
and a base point /o, the image of the monodromy representation 

M:7ri(P-{ai,...a n ,oo}) -> GL(2, C) 

if if /« 6\ (1-6) 

defined up to global conjugation by 

7 : (2/i,2/2)|/o = (2/i,2/2)|/o^7> ( L7 ) 
determines a subgroup C GL(2, C) that acts on the ratio 



Vi 

2/2 



(1.8) 



by linear fractional transformations (1.1). By a simple substitution of the type 

n 



i=l 



eq. (1.5) can be transformed into the form (1.3) without changing the projective class of the 
monodromy group; i.e., without changing the ©-action (1.1) on r. If (5 is Fuchsian and 
commensurable with SL(2,Z), the inverse function / = /(r) is a modular function in the 
above sense. 
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If R(f) is the resulting rational function in (1.3), then / = f(r) satisfies the Schwarzian 
differential equation 

{f,r} + 2R(f)f' 2 = 0, (1.10) 

where 

{/.->== t-1(t) 2 ' u ' :=d i ] (L11) 

denotes the Schwarzian derivative [HI, GS]. 

Perhaps the oldest example of such an equation involves the square of the elliptic modulus 
of the associated Jacobi elliptic functions A(r) = k 2 (r) : which satisfies the Schwarzian equation 
[HI] 

^ + ^w x ' 2 - - (L12) 

The automorphism group in this case is the level 2 principal congruence subgroup 

r(2) eSL(2,Z), = 1 (mod 2)}. (1.13) 
The associated Fuchsian equation is the hypergeometric equation 

, (1 _ A )^ + (1 _ 2A) |_l 1 , = 0, (1.14) 
for which a basis of solutions is given by the elliptic i-period integrals 



K -/V(i-«»Mi-^)-> i;1;t8) (1 ' 15a> 



iK' = J"-== == ^F(ii;l;l-0 (1.15b) 



eft 

v/(i - t 2 ~)(l - Pt 2 ) 2 

with 

r=— . (1.16) 

An equivalent way of representing the Schwarzian equation (1.12), due to Brioschi [B], is 
to introduce the functions 

1 d X' Id X Id, A' _ 

Wi:= 2Tr ln y W2:= 2^ ln (A^T)' " 3:= 2^ ln ACX^T)- (L17) 

These satisfy the system 

w[ = wi(w 2 + w 3 ) - w 2 w 3 

w' 2 = w 2 (w 1 + w 3 ) - wiw 3 (1-18) 

U>3 = Ws(wi + W 2 ) - W\W 2 i 
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introduced by Darboux [Da] in his study of orthogonal coordinate systems and solved by 
Halphen [Ha] with the help of hypergeometric functions. This system is referred to in [T, 01, 
02] as the Halphen equations. The general solution [Ha] is obtained by applying an arbitrary 
SL(2, C) transformation (1.1) to the independent variable, while transforming (2wi, 2w2, 2ws) 
as affine connections 

™<— > , ]_m w * oT C —7r ( L19 ) 

(cr + a) 2 ct + a 

The system (1.18) has appeared in a number of recent contexts, including: solutions of re- 
duced self-dual Einstein [GP] and Yang-Mills [CAC, T] equations, the 2-monopole dynamical 
equations [AH] , and the WDVV equations of topological field theory [Du] . 

A symmetrized version may be obtained by considering the symmetric invariants formed 
from (wi,W2,wz), such as 

W := 2( Wl + w 2 +w 3 ), (1.20) 
which satisfies the Chazy equation [C1,C2] 

W'" = 2WW" - 3W' 2 . (1.21) 

The corresponding modular function is obtained by noting that, although the action of the 
full modular group V := PSX(2,Z) upon A does not leave it invariant, the quotient T/r(2) 
by the normal subgroup T(2) is just the symmetric group S3, acting on A as the "group of 
anharmonic ratios" [Hi2] 

A^I.l-A,^,^,^!. (L22) 

These transformations correspond, respectively, to the following modular transformations 

r 1 1 t — 1 
r, ——, — , --— — , r+1, . 1.23 

1 + T T 1 + T T 

Symmetrization amounts to forming the ring of invariants, which has a single generator that 
may be taken as Klein's J-function 

j. Jfr'-A+l)' (124) 

whose automorphism group is Y. The corresponding Fuchsian equation is the hypergeometric 
equation 
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The associated Schwarzian equation, 



36 J 2 -41J + 32 



< J ' T > + J = °' (L26) 

was first studied by Dedekind [De] in relation to the modular properties of J and A. The 
resulting solution of the Chazy equation (1.21) is given [CI, T] by 

1 d ^ J' 6 d , A' 3 1 d , A 

where A is the modular discriminant. 

If we introduce analogous Halphen variables for the J-function 

^=^4 ^=^>(T^rr wr ' ! =5s to 7(^iy (1 ' 28) 

these satisfy the system 

1 1 23 31 31 

w{ =-w? + -wl + ^w x w 2 + ^w 1 w 3 - -W 2 W 3 

Wz 1 =\w* + X -Wl + ^W X W 2 - f^W 1 W 3 + ^W 2 W 3 (1.29) 

^3 = Vi 2 + \wi - |ViW 2 + ^W 1 W 3 + ^-W 2 W 3 . 
4 9 36 36 36 

From the transformations (1.22), it follows that the symmetrizing group S3 acts upon the 
Halphen variables (wi, w 2 ,w 3 ) by permutations. Differentiating the identity (1.24), we deduce 
that the quantities (Wi, W 2 , W 3 ) are related to (w\, w 2 , w 3 ) as follows 



3W 1 + 2W 2 + W 3 = 2(7! = W (1.30a) 

4(a 2 - 3a 2 ) 2 
2a f - 9a t a 2 + 27 a 3 



W 1 -W 3 = - n 3 n_ _ (l-30b) 



w w 2<T i - 9(71(72 + 27(7 3 n qo ^ 

W 2 -W 3 = 2 — » (1.30c) 

(T 1 — 3(T 2 

where 

o"! := + w 2 + W3, C2 := w\w 2 + w 2 w 3 + w 3 wi, a 3 := w\w 2 w 3 (1-31) 
are the elementary symmetric invariants. In terms of these, the system (1.29) reduces to 

(j'x = o 2 

a' 2 = 6a 3 (1.32) 
a- = 4aicr 3 - cr 2 
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(cf. [01]), which is equivalent to the Chazy equation (1.21). 

The modular solutions to these systems may also be represented as logarithmic derivatives 
of the null theta functions $2(7"), $3(1"), #4(7"), (cf. Appendix), in terms of which the elliptic 
modular function A(r) has the well known representation [WW] 



A(r) = 5 J = 1-3. (1.33) 

From the differential identities satisfied by $1, $3 (cf. Appendix, eqs. (A.19)-(A.21) and 
[01]), it follows that 

mt#2, v = i7n? i T7^-TT = Z7n? i ( L34 ) 



1-A 2 ' A 4 ' A(l-A) 

so the solution of the Halphen system may be expressed as 

wt = 2-^- In #4, w 2 = 2-^-lni? 2 , w 3 = 2-^-ln^3- (1-35) 
ar ar ar 

Using (1.24), the J-function may also be expressed rationally in terms of the theta func- 
tions 

(tf| + 0§ + 0|) 3 

Taking derivatives and applying the same differential identities, the J Halphen variables 
{Wi, W 2 , Ws) may be expressed as logarithmic derivatives of rational expressions in the #'s, 
and the symmetrizing relations (1.30a)-(1.30c) interpreted as differential relations satisfied by 
the theta functions. 

The above examples have been generalized by Ohyama [02] to other classes of second 
order equations of type (1.3), both Fuchsian and non-Fuchsian. For the Fuchsian case, the 
rational functions R(f) may be expressed as 

W = (1-3T) 

where 

n 

^(/)=n(/-°*) (i-38) 

i=l 

and N(f) is a polynomial of degree < 2n — 2. The generalization of the Halphen variables is 
given by 

X„:=iA ln/ , Xj:= Li ln _Z_ i= !,...„. (1 . 39 ) 
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For n > 2, the X^s again satisfy a system of first order equations analogous to (1.18) and 
(1.29), with suitably defined quadratic forms on the right, but they are also subject to a set 
of n — 2 quadratic constraints fixing the anharmonic ratios between any distinct set of four 
of them. The underlying phase space is therefore still 3c^dimensional and may, for generic 
initial conditions, be identified with the SX(2, C) group manifold. Besides the above two 
cases, other explicit solutions were given in terms of modular functions by Ohyama [03] for 
a system in which the corresponding Fuchsian operator has four regular singular points and 
the automorphism group is T(3). 

For the case of Fuchsian operators with three regular singular points the associated differ- 
ential systems were already studied by Halphen [Ha] . For the general hypergeometric equation 

/(I - /)^f + (c - (a + b + 1)/)^ - aby = 0, (1.40) 

the corresponding rational function is 

1/1-A 2 1-fi 2 \ 2 + ix 2 -v 2 -l\ , A , 

where 

A := 1 — c, n := c — a — 6, v := b — a (1-42) 

are the exponents at the regular singular points (0, l,oo), which determine the angles 
(A7r, /i7i~, utv) at the vertices. Introducing the variables 

1 d f Id f Id f 

w ^\tA- w ^\^jtt y Wtl= 2^-nrry (L43) 



and viewing them as functions of the ratio 



m 

V2 



of two independent solutions of (1.40), these satisfy the general Halphen system 

W[ =fi 2 W 2 + \ 2 Wl + v 2 Wl 

+ (z/ 2 - A 2 - [i 2 + l)WiW 2 + (A 2 - fi 2 - v 2 + l)WiW 3 + (/u 2 - A 2 - v 2 - 1)W 2 W 3 
W' 2 =\x 2 Wl + \ 2 Wl + v 2 Wl 

+ (u 2 - A 2 - fx 2 + l)WiW 2 + (A 2 -[i 2 -v 2 - l)WiW 3 + (fi 2 - A 2 - v 2 + 1)W 2 W 3 

Wz =[i 2 Wl + \ 2 Wl + v 2 Wl 

+ (u 2 - A 2 - fx 2 - 1)W!W 2 + (A 2 - \i 2 - v 2 + l)WiW 3 + (/u 2 - A 2 - v 2 + 1)W 2 W 3 . 

(1.44) 
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Although this gives a construction, in principle, of the general solution to such systems, 
the functional inversion involved is not in general globally well defined. Only if the parameters 
(A, fj,, u) are of the form ( — , -, |), where m, n and p are integers or oo, is there a tessellation 
of the upper half plane by (triangular) fundamental regions, and even then we cannot in 
general say much about the explicit form of the inverse function. However, for a small number 
of special cases which are described below, these again turn out to be modular functions 
that can be given explicit rational expressions in terms of null ^-functions or the Dedekind 
^-function. 

Remarkably, the J-function also appears in connection with "Monstrous Moonshine" 
[CN], in that the q = e 27riT -coefficients of the Fourier expansion of j := 12 3 J — 744 are the 
dimensions of representations of the Monster simple sporadic group; i.e., the traces on the 
identity element. Like j, the (/-series with coefficients given by traces on the other conjugacy 
classes turn out to also be Hauptmoduls; i.e., each is the generator of a field of meromorphic 
functions of genus 0. An alternative characterization of the elliptic modular function j is the 
fact that the principal part of its g-series is g _1 , together with its behaviour under the action 
of the Hecke operator; namely, 



where P n j(f) is the Faber polynomial [Fa] of degree n. The functions appearing as such 
character generators are included in a larger class of Hauptmoduls, the replicable functions 
[CN, FMN] , which are constructed using a generalization of the Hecke operator. Their auto- 
morphism groups all contain a finite index subgroup of the type 



(and hence their automorphism groups are all commensurable with V.) The maximal such N 
is referred to as the level of the function. Like all Hauptmoduls, they also satisfy Schwarzian 
equations of the type (1.10) with rational R(f) of the form (1.37), and hence each has an 
associated Fuchsian operator of the form (1.3) whose projectivized monodromy group is the 
automorphism group of the function. 

In the present work, we study the differential equations satisfied by such Hauptmoduls and 
their corresponding generalized Halphen systems. In particular, we consider the Hauptmoduls 
known as triangular functions, for which the tessellation of the upper half-plane is associated 
with a triangular domain and a corresponding hyperbolic group generated by reflections in 
its sides. The associated Fuchsian equations are therefore of hypergeometric type. Since the 
automorphism groups appearing are conjugate to subgroups of PGL(2, Q), it follows that only 




Vn > 1, 



(1.45) 



ad — rt 




(1.46) 
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crystallographic angles (0, 7r, 7r/2, 7r/3, 7r/4, 7r/6) occur at the vertices. For the j-function, the 
triangle has a cusp at zoo and angles of 7r/3 at e l7r / 3 = J _1 (0) and tt/2 at i = J _1 (l), so this is 
denoted (|, ^ 0). Since the functions considered have (/-expansions with principal part (/ _1 at 
q = 0, there must be a cusp at r = zoo, and hence the angle v always vanishes. There are nine 
such triangular arithmetic groups [Ta] with angles given by: (0,0,0), (1/2,0,0), (1/3,0,0), 
(1/3, 1/2, 0), (1/4, 1/2, 0), (1/6, 1/2, 0), (1/3, 1/3, 0), (1/4, 1/4, 0), (1/6, 1/6, 0), and these are 
precisely the ones appearing in the list [CN] of replicable functions with integral (/-series 
coefficients. 

In the following sections, several examples will be given of differential systems of general- 
ized Halphen type whose solutions are completely determined in terms of replicable functions. 
These will include all the arithmetic triangular groups, and some further cases for which there 
are four singular points in the associated Fuchsian equation (as well as one in which there are 
26). Each such function has a fundamental domain bounded by circular arcs centered on the 
real axis, with a cusp at ioo and one of the crystallographic angles (0, n/2, n/3, 7r/4, ir/6) at 
each remaining vertex. By construction, each has a normalized (/-series of the form [FMN] 

with integer coefficients a n . In the case of triangular functions, in order to relate these 
to solutions of the hypergeometric equations, another standard normalization is chosen by 
applying an afhne transformation 

F = af + p, (1.48) 

with constants a and (3 chosen to assign the values (0, 1, oo) to f(r) at the vertices of the 
fundamental triangle. Since all these functions are known explicitly [FMN] in terms of the 
Dedekind eta function rj(r) or the null theta functions ^(t), we are able to give explicit 
solutions of the corresponding differential systems in terms of logarithmic derivatives of the 
rj- or ^-functions. 

Certain of these functions are expressible in terms of others through a rational map, 
possibly composed with a transformation of the type r i— > r jn (where n has prime factors 2 
or 3). This implies that there is a subgroup relation between their respective automorphism 
groups, and also that the differential system associated to the larger automorphism group 
is a symmetrization of the one associated to the subgroup. An example of this was seen 
in the pair of elliptic modular functions (J, A), where the rational expression (1.24) implies 
the subgroup relation T(2) C V and the relations (1.30a)-(1.30c) between the correponding 
Halphen variables. In general, if / satisfies the Schwarzian equation (1.10) and g satisfies 

{g,T} + 2R(g)g' 2 = 0, (1.49) 
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then if / can be expressed as a function of g 

f = Q(s), (1-50) 

this will satisfy the Schwarzian equation 

{Q,g} + 2R(Q(g))Q' 2 = 2R(g). (1.51) 

Equivalently, if Q(g) satisfies (1.51) and y(f) is a solution to (1.3), the composite function 
y(g) = (Q')~h(Q(g)) will satisfy 

fl+R(g)y = 0. (1.52) 

By choosing both R(f) and R(g) of the hypergeometric form (1.41), Goursat [Go] found ra- 
tional transformations of degree < 4 relating various classes of hypergeometric functions. The 
relation (1.24) may be viewed as a composition of two such transformations, of degrees 2 
and 3, connecting the hypergeometric equations of types (a, 6, c) = (j^, j^, |) and (|, |, 1). 
In section 2b, transformations of Goursat's type are used to relate the associated arithmetic 
triangular functions, and in section 3b, similar transformations of degree < 4 are given relat- 
ing these to certain 4-vertex cases. Such relations are determined by finding appropriately 
normalized rational solutions / = Q(g) of the Schwarzian equation (1.51) or, equivalently, by 
equating the (/-series of both sides up to a finite number of terms. 

In all the cases treated, we give a tabular summary of the properties of the modular 
functions /, including explicit rational expressions for them in terms of the rj- and ^-functions. 
The general solutions to the Halphen systems (1.44) and their multivariable generalizations are 
thus determined as logarithmic derivatives of rj- or ^-functions. For the hypergeometric cases 
(section 2a), we list the exponents (A, /i, v) determining the Halphen type system (1.44), the 
corrresponding hypergeometric parameters (a, 6, c) and the group elements fixing the vertices. 
For the four vertex cases (section 3a), we give the locations (a±, a 2 , a 3 ) of the finite poles of the 
rational function R(f); i.e., the values of / at the vertices, the group elements stabilizing the 
vertices and a representative of an associated linear class of quadratic forms in the dynamical 
variables that serves to uniquely define the corresponding constrained 4-variable system, as 
well as the rational function R(f) entering eqs. (1.3) and (1.10). Section 4 contains a discussion 
of the n + 1 pole case and includes an example with 26 regular singular points, which is the 
largest number that appears, and automorphism group of level 72. In the general n + 1 pole 
case, it is shown that Ohyama's quadratically constrained n + 1 variable dynamical system 
[02] is equivalent to a 3-variable system defined on the SL(2, C) group manifold. 

For both the triangular cases (section 2b), and certain four vertex cases (section 3b), we 
examine the various rational maps of degree < 4 relating the different Hauptmoduls, together 
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with the corresponding relations between the associated automorphism groups. These maps 
define symmetrizations for the function fields and the corresponding differential systems. Just 
as in the J- A case (1.30a)-( 1.30c), the generalized Halphen variables for the more symmetrical 
systems consist of rational invariants formed from the less symmetrical one. Not all cases 
correspond to normal subgroups, however, so the fibre of the quotient map is not necessarily a 
Galois group, but rather the quotient of two Galois groups, corresponding to a field extension 
whose automorphism group is the largest subgroup normal in both the groups of the pair of 
Hauptmoduls. 

2. Triangular replicable functions and symmetrization maps 

2a. Solutions of Halphen type systems. 

In Table 1 , we list all cases where the underlying Fuchsian equation is of hypergeometric 
type (1.40). The Halphen system is always of the form (1.44), and therefore we list only the 
values of the constants (a, b, c) and (A, p,, v) characterizing these systems. The notation used 
to designate the corresponding automorphism groups and function fields is that of [FMN]. 
The triangular cases appearing are denoted: 1A, 2A, 3A, 2B, 3B, 4C, 2a, 4a and 6a, where 
the integer N denotes the level. (The upper case letters denote functions that are character 
generators for the Monster.) There are further cases of replicable functions having triangular 
automorphism groups, but these are all equivalent to one of the above under an affine transfor- 
mation that just relocates the two finite poles of the Fuchsian equation (1.3), together with a 
Mobius transformation (1.1) in r, and therefore they are not listed separately. The ^-function 
formulae listed for the functions /(r) normalized to take values (0, 1, oo) at the vertices are 
obtained from those given in [FMN] for the functions F normalized as in (1.47) by apply- 
ing the appropriate affine transformation (1.48). The ^-function formulae are deduced either 
from standard identities relating the rj- and functions (cf. Appendix), or by applying the 
rational maps relating the different cases listed in sections 2(b) and 3(b). The automorphism 
group, being the projective image of the monodromy group of the associated hypergeometric 
equation, has three matrix generators po, p\, p^ stabilizing the vertices satisfying the relation 



for some m ^ 0. Representatives of the projective class of these generators may always be 
chosen to have integer entries, though not necessarily unit determinant, so the constant m in 
(2.1) is an integer. In all cases the generator p^ stabilizing the cusp at zoo is 



PooPiPo = ml, 



m G R 



(2.1) 




(2.2) 



and the two remaining elements po, p\ are given in the table. 
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Table 1. Triangular Replicable Functions 



Name 


(a, 6, c) 


(A, fi, v) 


Po Pi 


rp 

r 


J{ T ) 


1A 

~ r 


(± JL ^ 

V 12' 12 ' 3^ 




c=o (?-o 


1728/ - 744 




2A 


v 8 ' 8' 4^ 




(2-2) (2 J) 


256/ - 104 




3A 


V6' 6' G> 




(3 -3) (° 0) 


108/ - 42 


(r 7 12 (r)+27 ?7 12 (3r)) 2 
108t7 12 (t)?7 12 (3t) 


2B 

~r (2) 


(1 I I) 

V4' 4' 2/ 




A -A /-1 o\ 
{2 -1) \ 2 -1) 


64 f - 40 


1,1/ »7(t) ^i 24 
1 "r" 64 Vr 7 (2r)>' 

( 1 ? 4 (r)+<(r)) 

0»(t) 


3B 

~r (3) 


C- - -) 

V3' 3' 3' 


(|,0,0) 


( 1 —i\ ( — 1 o N \ 
I 3 -v \ 3 -v 


27/ - 15 


1 , 1 / V(r) \ 12 
L ^ 27\r)(3T)J 


AC* 

~r (4) 


( l ~ 1) 

\ 2 ' 2 ' / 


(0 0) 


A -A (-1 0) 


16 f - 8 


1 *S(2t) 
A(2r) - 0*(2t) 

- 1+ 1 C ^ I s 

T 16 Vt?(4t)/ 


2a 


V6' 6' 3^ 


V 3 ' 3' u / 


/ 2 _ 3 \ /o -A 

I 4 -v I 4 -v 


24v / 3zf2f - 1) 


V / 3 l (e^/ 3 1 ? 4 (2r)- 1 ? 4 (2r)) 3 
9^(2r)^(2r)^(2r) 


4a 


(I I 3n 

\ 4 ' 4' 4^ 


(5,5,o) 


(s _s) (s 4) 


-16z(2/-l) 


l (^ 2 (2r)+i 1 9l(2r)) 4 
8 1 9 4 (2r) 1 ? 2 (2r)^ 2 (2r) 


6a 


V3' 3' 6^ 




\^12 -12 y \12 -6) 


6^(2/ - 1) 


V / 3i(?7 6 (2r)+3V3i?7 6 (6r)) 2 


36?7 6 (2t)7/ 6 (6t) 



* Remark: Case 4C has the same fundamental domain as the elliptic modular function A(r), but with re- 
spect to the variable t/2. The automorphism group is r (4), which is conjugate to r(2) under the map 



2b. Rational maps between triangular cases. 

In the following, we catalogue the rational maps, analogous to (1.24), relating the repli- 
cable functions listed above. Each map may be seen as following from an identity relating 
different hypergeometric functions through a rational change in the independent variable. Such 
identities were derived systematically by Goursat [Go], and are labelled here according to the 
numbering scheme of [Go]. For each case, we have a pair (/, g) of function field generators, 
with / a rational function of g, normalized consistently with the hypergeometric equations, 
so that the vertices of the fundamental triangle are mapped to (0, 1, 00). In order to preserve 
the standard ordering with angles (A, fj,, v) at the vertices with values (0, 1, 00), we make use 
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of the fact that F(a, 6, c; z) and F(a, 6, 1 — c + a + b; 1 — z) satisfy the same equation and 
compose, in some cases, with the maps / i— > 1 — f or g i— > 1 — g. The corresponding pair of 
functions, normalized as in (1.47), is denoted (F, G) and the relevant map is denoted G i— > F. 
In order to preserve the form (1.47) of the g-series, it is necessary in some cases to compose 
the rational map with a map r i— > ^ or r i— > ^ on the independent variable. In these cases, the 
composed map is denoted C i— > F. By taking derivatives of the corresponding rational map 
relating / to g, we also obtain relations between the Halphen variables, such as that given in 
the case of the pair (J, A) by (1.30a)-(1.30c). 

To fix notation, the Halphen variables for the two cases will be denoted 

1 d f Id f Id f 

W »=2**T w '-=2^jrry w ^2>nh (2 ' 3) 

for the first element of the pair (/, g), and 

1 d , g' Id, s' 1 d , g' . 

for the second. Since / is in the function field generated by g, it is (J5 ff -invariant; i.e., <3 g C 
<5f. When <5 g is a normal subgroup, <3f acts on g by rational tranformations, and on the 
Halphen variables (w\, W2, W3) linearly. Since this defines a finite group action, the ring of 
invariants in the (tui, W2, W3) variables is generated by 3 elements, and the (Wi, W 2 , W3) are 
uniquely determined as rational expressions of these. In all the cases listed below we give 
the hypergeometric identity underlying the rational maps between the modular functions and 
the maps, expressed in terms of both the (f,g) normalizations and the (F,G) ones. When 
C5 S C <3f is a normal subgroup, we give the symmetrizing group := <3f/<3 g (i.e., the 
Galois group of the field extension), the action of (3f upon g (analogous to the S3 action 
(1.22) on the elliptic modular function A(r)), the elementary symmetric invariants under this 
action, and the linear relations determining (Wi, W2, W3) in terms of the latter following from 
differentiation of the rational map relating / to g. When (& g C 0/ is not normal, we list a 
symmetrization quotient Sj/S^, where and are the symmetrizing groups of the smallest 
field extension, with generator h, that is Galois for both. The symmetrizing groups all turn 
out to be either a symmetric group S n , an alternating group A n or a cyclic group Z n . 

We note that in each case, there exists at least one relation involving a linear symmetric 
invariant of the form 

PW 1 + QW 2 + RW 3 = pw 1 + qw 2 + rw 3 (2.5) 
where (F, Q, F,p, q, r) are all integers, and 

P + Q + R = p + q + r=:k. (2.6) 



13 



J. Harnad and J. McKay 



This is deduced from a relation of the form 



j'P+Q+R g'p+q+r 

(2.7) 



f p +R(f - i)Q+R g p+r (g - l) q+r ' 

where (P, Q, R,p, q, r) are integers, which may always be found by simply choosing (P, Q, R) 
so that the left hand side of (2.7) has no singularities at the vertices. The resulting quantity 
is an analytic form of weight 2k, and its logarithmic derivative, given by (2.5), satisfies a third 
order equation analogous to the Chazy equation (1.21). 

2a' i-> 1A: 

Hypergeometric identity (Goursat (138)): 

112. ,\ (\ 1 2 \ . . 

— ,— :-:4x(l-x) )=F i-,-:-:x) (2.8) 
12' 12' 3' V ') U'6'3' j V ; 



Rational map: 

F = 984 + G 2 o (r h-> r/2), 
/ = 4</(l - </) o (r -> t/2) 
Symmetrizing group: S 2 ^ = Z 2 



(2.9) 



r g(r/2)^ 1 - g(r/2), (w 1 ,w 2 ,w 3 ) ^ (w 2 ,w 1 ,w 3 ) (2.10) 

r 

Polynomial invariants: 

a 1 := Wi + w 2 , a 2 :=wiw 2 , Ei := ty 3 (2.11) 
Powers and coefficient in (2.7): 

(P, Q, R- p, q, r; M) = (3, 2, 1; 2, 2, 2; -1) (2.12) 
Relation between the Halphen variables: 

Wi + W 3 = Si o (r i-> r/2) 

W\ - W3 = — _ — o (r ^ r/2) ( 213 ) 

W 2 - W 3 = ]-a x - Ei o ( T m r/2) 

4a' i-> 2A: 

Hypergeometric identity (Goursat (138)): 
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Rational map: 

F = 152 + G 2 o(r i->t/2), 
/ = 4</(l - </) o (t -> t/2) 

Symmetrizing group: = Z2: 

r — , 5r(r/2) l-#(r/2), (^1,^2,^3)^(^2,^1,^3) 
At 

Polynomial invariants: 

(Ti := w\ + W2-, 02 := W1W2, Ei := w 3 
Powers and coefficient in (2.7): 

(P,Q,R;p,q,r;M)=(2, 1,1; 1,1,2;-^) 

Relation between the Halphen variables: 

Wi + W 3 = £1 o (r i-> r/2) 

Wi - W3 = — o (r ^ r/2) 

O"! — ZZji 

W 2 - W 3 = -o-i - Si o ( T m r/2) 

6a' i-> 3A: 

Hypergeometric identity (Goursat (138)): 
Rational map: 

F = 66 + G 2 o (r i-> r/2), 
/ = 40(1 -(/)o(t- r/2) 
Symmetrizing group: S^' = Z 2 : 

r^-^-, o(r/2) l-o(r/2), (w 1 ,w 2 ,w 3 )^(w 2 ,w 1 ,w 3 ) 
Polynomial invariants: 

(7 1 := W\ -\- W2, (72 := W1W2, Ei := ws 
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Powers and coefficient in (2.7): 



(P,g, J R;p,g,r;M)=(3,l,2;l,l,4;^) 

lb 



Relation between the Halphen variables: 

Wi + W 3 = Ei o (r h-> r/2) 



- W3 = — o (r ^ r/2) 

G\ — ZZji 

- W 3 = -o-i - Ei o (r i-> r/2) 



4C i-> 2S: 

Hypergeometric identity (Goursat (44)): 



Rational map: 



,1 1 x 2 \ (\ 1 

(] - :rTJ -1.4'J :li,:= i(^) =F U'2 !l!I 



f = c + 7 ^. /- (9 + 1)2 



G + 8' ' ~ Ag 
Symmetrizing group: S^b = ro(2)/r (4) = Z 2 : 

T 1 > _ T , t , </•-»■ -, (tUl, ^2,^3) >-> (Wl,^,^) 

2r + 1 g 

Polynomial invariants: 

0"i := w 2 + w 3 , cr 2 := W2W3, Ei := w\ 
Powers and coefficient in (2.7): 

(P, Q, R; p, q, r; M) = (2, 1, -1; 4, -1, -1; 1) 
Relation between the Halphen variables: 

W 1 + W 2 = 2Ei 
Wx - W 3 = 2Ei - ai 
<j\ - 4a 2 



W 2 - W 3 = 



2Ei - (7i 
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AC h-> 2B: 

Hypergeometric identity (Goursat (138)): 



F Q, 1; l;4r(l - *)) = F Q, i; l;z) (2.32) 



Rational map: 

F= -40 + G 2 o(r i->r/2), 
/ = (2<7 - l) 2 o (r ^ r/2) 

Symmetrizing group: S|g' = r (2)/r ( 2) = Z 2 : 



(2.33) 



n->r+l, </(r/2) i-> l-g(r/2), (w 1 , w 2 , w 3 ) i-> (w 2 , W3) (2.34) 
Polynomial invariants: 

(7i := toi + mj 2 , a 2 \=WiW 2 , £1 := iy 3 (2.35) 
Powers and coefficient in (2.7): 

(P, Q, fl; p, g, r; M) = (2, 1, -1; 2, 2, -2; 4) (2.36) 
Relation between the Halphen variables: 

W 2 + W 3 = Si o (r i-> r/2) 

Wi - W3 = ^1 -Eio(th r/2) (237) 
lT a -i y , = ^- E '°i +S ') .(r»r/2) 

2B ^ 2A: 

Hypergeometric identity (Goursat (44)): 

,i/113 x 2 \ /111 \ , 
(1 - z)sF -, -; -; z := = F -, -; -;x (2.38) 



8'8'4' ' 4(x-l)/ V 4 '4'2 
Rational map: 

2 12 o 2 
F = G + CT '=4&3Tj < 2 ' 39 > 

Symmetrizing group: S 2 f = r[, 2) (2)/ro(2) = Z 2 : 

r — , (/!-»• — 9 — , (w 1 ,w 2 ,w 3 ) ^ (w 3 ,w 2 ,w 1 ) (2.40) 
2t g - 1 
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Polynomial invariants: 

a\ := w\ + w 3 , a 2 :=wxw 3 , Ei := w 2 (2-41) 

Powers and coefficient in (2.7): 

(P, Q, R- p, q, r; M) = (2, 1, 1; 1, 2, 1; 4) (2.42) 

Relation between the Halphen variables: 

Wx + W 2 = 2Ei 

Wl -W 3 = (2.43) 
W 2 -W 3 = -a x + 2Ei 

3B ^ 3A: 

Hypergeometric identity (Goursat (44)): 

(1 - -, -; -;z := — r = F ( -,-;-; a; ] 2.44 

V ; V 6 6 6 4(x-l)J \3 '3'3' y V ; 

Rational map: 

Symmetrizing group: Sg^ = Z 2 : 

r — , — ^-r, (twi,tW2, w 3 ) ^ (^3,^2,^1) (2.46) 
3r # - 1 

Polynomial invariants: 

(Ti := Wi + w 3 , a 2 :=w l w 3n Ei := w 2 (2.47) 
Powers and coefficient in (2.7): 

(P, Q, R- p, q, r; M) = (3,1, 2; 2, 2, 2; 16) (2.48) 
Relation between the Halphen variables: 

W x + W 2 = 2Ei 

W 1 -W 3 = (2-49) 
- W 3 = -ax + 2Ex 
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IB i — > 1A: 

Hypergeometric identity (Goursat (123)): 

„ / 1 11 x(x - 9) 2 \ , , i (1 1 1 

^ — , — ;-; ^ = (1 ( -, -; -;x 

V12' 12' 2' 27(x-l) 2 / V ; V 4 4 2 

Rational map: 

2 16 (3G + 184) 



F = G + 
/ = 



(G - 24) 2 
(<7 + 3) 3 



270/ -1)2 

Symmetrization quotient: T/ro(2) = Sf^/S|^ = S3/Z2 
Polynomial invariants: 

Ei := Swi + 2w 2 + W3, E 2 := (wi - w 3 )(4wi - 3w 2 - 
E 3 := (wi - ws) 2 (8wi - 9w 2 + W3) 

Powers and coefficient in (2.7): 

(P, Q, R; p, q, r; M) = (3, 2, 1; 3, 2, 1; -27) 
Relation between the Halphen variables: 

2Wi + 3W 2 + W 3 = Ei 

E 2 

^-^ = --2 

^3 

E 3 

2S' 1 — > 1A: 

Hypergeometric identity (Goursat (122)): 

Vl2'l2'2' 27(1 -x) J [ ] V4' 4 '2 

Rational map: 

F = G 2_ 552 + _2^_ o(r ^ r/2)5 



27(1-^) 1 1 ' 



Symmetrization quotient: T/ro(2) = Sf^ /S 2 b = S3/Z2 
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Polynomial invariants: 

Ei := 3wi + w 2 + 2w 3 , E 2 := (wi - w 3 )(wi + 3w 2 - 4w 3 ), 
E 3 := (wi - w 3 ) 2 (wi - 9w 2 + 8w 3 ) 

Powers and coefficient in (2.7): 

27 

(P, Q, R; p, q, r; M) = (3, 2, 1; 3, 1, 2; - — ) 



64' 



Relation between the Halphen variables: 



3Wi + 2W 2 + W 3 = -Ei o (r ^ r/2) 

W 1 -W 3 = -^o(r^r/2) 
Zl^ 3 

W 2 -W 3 = -^o( t ^t/2) 

ZXjo 



3B ^ 1A: 

Hypergeometric identity (Goursat (131)): 



,1 1 2 x(x + 8) 3 \ ,„ ,i (\ 1 2 
( To' To' o S cTTT TTs = (1 - ^ o' o 5 ^ * 



12' 12' 3' 64(x- 1)V v ' V 3 ' 3 ' 3 
Rational map: 

3 9 (10G 2 + 732G + 9459) <7(<7 + 8) 3 

^ = G H 77^ 77^ ' J 



(G-12) 3 ' J 64(^7- l) 3 

Symmetrization quotient: r/r (3) = Sf|/S§| = A 4 /Z 3 
Polynomial invariants: 

Ei := 3w\ + 2w 2 + w 3 , E 2 := {w\ - w 3 )(9w 1 - 8w 2 - w 3 ), 
E 3 := (wi — w 3 ){27w\ — 36w\w 2 + 8w 2 — 18wiw 3 + 2w 2 w 3 — w 3 ) 

Powers and coefficient in (2.7): 

{P,Q,R;p,q,r;M)= (3, 2,1; 3, 2,1; 64) 
Relation between the Halphen variables: 

3Wi + 2W 2 + W 3 = Ei 

E 2 

^-^ = --2 

^3 

s 3 

^2-W 3 = --^ 
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3B' i — * 1A: 

Hypergeometric identity (Goursat (130)): 



1 1 2 x(9x-8) 3 \ ,i/112. 

1 (2.65) 



12' 12' 3' 64(x- 1) 7 v y V 3 ' 3 ' 3 



Rational map: 



F = G 3 - 162G + 228 + f - 744 o (r r/3), / = ^ o (r r/3) (2.66) 

Symmetrization quotient: r/r (3) = S^'/ S ll' = A 4/Z 3 
Polynomial invariants: 



Ei := wi + 2w 2 + 3w 3 , E 2 := (wi - w 3 )(w 1 + 8w 2 - 9w 3 ), 
E 3 := (wi — ws)(wf — 20wiW2 — 8W2 + I8W1W3 + 36W2W3 — 27m; 2 ) 



Powers and coefficient in (2.7): 



(2.67) 



64 

(P, Q, R- p, q, r; M) = (3, 2, 1; 2, 3, 1; — ) (2.68) 



Relation between the Halphen variables: 

3Wi + 2W 2 + W3 = ^Ei o (r ^ r/3) 

E 2 

- W 3 = --f- o (r h+ r/3) (2.69) 
W 2 -W 3 = -^r-o(r^r/3) 

AC i-> 2a: 

Hypergeometric identity: 

F (-,-;-; — 7= — ) =(x(l-x))*F [ -,-;l;x), u := e~ 2.70 

V6' 6' 3' sVsx(i-x)J 1 1 " V 2 2 1 ' 

Rational map: 

f - G+ G?36i- '- 3v ^(l- 9) (2 ' 71) 
Symmetrizing group: = Z3: 

t >-» ~ . n , (wi,«;2,W3) ^ (^3,^1,^2) (2.72) 

4r - 2 g 
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Polynomial invariants: 
Powers and coefficient in (2.7): 



= Wi + W 2 + Ws, 

= wiw 2 + w 2 w 3 + wiw 3 , (2.73) 



(P, Q, R; p, q, r; M) = (1, 1, 1; 1, 1, 1; 3v^) (2.74) 
Relation between the Halphen variables: 

W 1 +W 2 + W 3 = °1 

(a? - 3a 2 ) 2 



W 2 -W 3 = — 



S 3 (2-75) 
S 3 



o\ - 3(T 2 

3. Four vertex systems 

3a. Solutions of Generalized Halphen systems. 

The examples listed in Table 2 below all involve Fuchsian operators (1.3) with four regular 
singular points (including 00). We use the following linear combination of Ohyama's dynamical 
variables [02] as our phase space coordinates 

u:=X = ^, v r .= ^{X Q -X i )= 1 -j^— i = l,2,3, (3.1) 

where (01,02,03) are the locations of the finite poles of the rational function R(f) in the 
Fuchsian equation (1.3) and the Schwarzian equation (1.10). The cases considered here involve 
the replicable functions denoted: 6C, 6D, 6E, 6c, 6E, 9B in [FMN]. These include all cases 
with integer (/-series coefficients that, when composed with a suitably defined rational map 
of degree < 4, give one of the triangular functions listed in Table 1. Each provides modular 
solutions to a system of the type formulated by Ohyama [02], generalizing the equations of 
Halphen type. Further cases of such replicable functions do exist, but they may all be related 
to one of the above through an affine transformation that relocates the finite poles of the 
Fuchsian operator (1.3), composed with a Mobius transformation of the modular variable r, 
and hence they are not listed separately. As in the previous section, the function F(r) denotes 
the normalized g-series (1.47) as in [FMN], which is related by an affine transformation (1.48) 
to a function f(r) normalized, if possible, to take values (a, 0, 1, 00) at the vertices, for some 
real a < 0. In two cases: 6D and 9B, this is not possible, because the three finite singular 
points are not collinear. In case 9B, they are normalized instead to the cube roots of 1, as in 
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[03], while in case 6D they are normalized so that the single real pole is located at 1. We also 
list, for each case, expressions for /(r) in terms of the ry-function or null ^-functions and the 
generators {pi, p 2 , P3} of the automorphism group corresponding to the finite vertices of the 
fundamental region normalized to have integer entries. These satisfy 

Poopzp2p\ = ml (3.2) 

where is given by (2.2) and m 7^ is an integer. 

The differential equations satisfied by the quantities (u, v±, v 2 , ^3) are determined by a 
quadratic vector field together with a quadratic constraint. The equations for v[,v 2 ,v' 3 fol- 
lowing from (3.1) are of the form 



vl = -2vf + 2uv h i = l,2,3, (3.3) 
while the constraint equation is 

(ai - a 2 )viv 2 + (a 2 - a 3 )v 2 v 3 + (03 - ai)viv 3 = 0. (3.4) 
When the rational function R(f) appearing in eqs. (1.3) and (1.10) is expressible in the form 

L 1J — X 



the remaining equation, for u', is 



u' = u 2 



Yl r iJ v i v r ( 3 - 6 ) 



In the table, rather than listing R(f) directly, we give the quadratic form in (vi,V2,vs) 
appearing on the RHS of eq. (3.6). (Note that this quadratic form, and hence also the coeffi- 
cients Tij appearing in (3.5), is arbitrary up to the addition of any multiple of the constraint 
(3.4), but such a change leaves R(f) invariant.) The general solution is again obtained by 
composing the function / with a Mobius transformation (1.1). This amounts to transforming 
the function 2u as an affine connection 

u — ► -. 1 -—,uoT — v (3.7) 

(cT + d) 2 cr + d v ; 



and the functions Vi as 1-forms 



Vi^ , , ,^ i°T. (3.8) 
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Table 2. Four Vertex Replicable Functions 

(admitting a rational map of degree <4 to a triangular replicable function) 



Name 


(ai,a 2 ,a 3 ) 


Pi 

P2 
P3 


3 


F 


fir) 


6C 


(-3,0,1) 


I 12 -3 ; 

\ / 

(1 


) 


-|v 2 f3 - ^1^3 


4/ + 2 


i , 1 -) a (r)r) 6 (3r) 
1 "r 4 77 6 (2t)?7 6 (6t) 


6D 


GM,i) 

/3:=-| + v / 2^ 


(e J) 

(;=0 
(.-; 


) 


3^1 + ~A + vi 

28 -J?^^1^3 
28+ ^^3 




1 1 r, 4 (r)^(2r) 
1 "r" 4 ?7 4 (3t)?7 4 (6t) 


6E 

~ r (6) 


(-|,o,i) 


(i 5 2 -Vj 

( 5 -2> 
I 18 -7 i 

\ / 

(11 


) 


+ v\ + vl 
-^v 2 v 3 - %ViV 3 


8/ -3 


i | 1 7? 5 (r)77(3r) 
^8 ?7(2r)?7 5 (6r) 


6c 


(-1,1,0) 


(8-7 
\12 -10 

/ 2 -iS 
I 12 -4 ; 

\ / 

Ch zL 


) 
) 


+ lfv 1 v 2 




3V3 ry 6 (6T) 


8E 

~ r (8) 


(-1,0,1) 


^ 16 -5y 


) 


vf + v% + vi 

-2V! v 3 


4/ 


111 7 / ( T )V ) 
1 ~T 4 ?7 2 (2t)?7 4 (8t) 

0§(2t)+*2(2t) 
^(2r)-^(2r) 


9B 

~ r (9) 


(w,a>, 1) 

2-iri 

w :=e 3 


(i -i) 


) 


v I + v I + vi 

—V1V2 

-(1 - u)vxv 3 
-(1 - u)v 2 v 3 


3/ 


1 , i 

1 3 773(9r) 



3b. Rational maps to triangular cases. 

In the following, we catalogue the irreducible rational maps of degree < 4 which, like those 
in Section 2b, relate the modular functions listed above to the triangular cases of Section 2a. 
The same notation is used as in Section 2b, (F, G) denoting a pair of replicable functions, with 
F given as a rational function of G, where F is one of the triangular functions of Section 2a, 
with normalized (/-series given by (1.47), and the four vertex function G is one of the cases 
listed in Section 3a, similarly normalized. The corresponding pair, with values normalized to 
(0, 1, oo) at the vertices in the triangular case and to the values (ai, a 2 , a 3 ) listed in Table 2 
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for the four vertex case, is denoted (f,g). The generalized Halphen variables for / are again 
as defined in (2.3), while the corresponding ones for g are defined by replacing / in (3.1) 
by g. Taking derivatives of the rational map relating / to g, we obtain relations between 
the associated generalized Halphen variables. Analogously to relation (2.7), there is always a 
relation of the form 

fiP+Q+R ik 

= M- — - — — -, (3.9) 



fP+R(f - 1)Q+R ( g - ai )P(g - a 2 y{g - a 3 )' ' 
with integer powers (P, Q, R, p, q, r, k) satisfying 

P + Q + R = k 7 (3.10) 

and (.P, Q, R) again chosen so that the LHS of (3.9) has no singularities at the vertices. The 
resulting quantity is again an analytic form of weight 2k and, taking logarithmic derivatives of 
both sides of (3.9), we again obtain a linear relation between the generalized Halphen variables 
for the pair (/, g) 

PWi + QW 2 + RW 3 = ku- pvi - qv 2 - rv 3 . (3.11) 

Two other relations also follow from the definitions of these variables, allowing us to determine 
the quantities (W\, W 2 , W3) for the triangular case as simple rational symmetric functions of 
the variables (^,^1,^2,^3) for the four vertex case, invariant under the larger automorphism 
group (3f. In those cases where <3 g C (25/ is a normal subgroup we list the linear, quadratic 
and cubic polynomial invariants, in terms of which (Wi, W 2l W3) are expressed. In all cases a 
fourth invariant, which we do not list separately, is provided by the vanishing quadratic form 
(3.4). (A notational convention that is slightly different from that of Table 2 is used below; 
the subscripts on the f-variable designate the location of the poles rather than their order in 
the sequence (1,2,3) indicated in the second column of the table; i.e., what appears in the 
table as vi is here denoted v ai .) 



8E 1 — > AC: 
Rational map: 



2 4 

F = G+ G' 



(3.12) 



4S 

Powers and coefficient in (3.9): 

(P, Q, R; p, q, r, k; M) = (1, 0, 0; -1, 1, 1, 1; 1) (3.13) 
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Symmetrizing group: S|g = r (8)/r (4) = Z 2 : 

-r 1 

T I— > , OH -, 

4r-l' y a' (3.14) 
(it, v_i,v , vi) h-> (u-2v ,v 1 -v ,-v ,v 1 -v ) 

Polynomial invariants: 

Ei := u + v-i — v — vi, E' x := 4ui — 2i? , E 2 := 4i>q (3.15) 
Relation between the Halphen variables: 



W t = Ei 
Ej 
E' 



- W 3 = |£ . (3.16) 



W 2 -W 3 = E'i 



6c i — > 6a: 
Rational map: 



3 3 

f= (<7 + l) 2 
J " 4<7 

Powers and coefficient in (3.9): 



Wi = Ei 
W 1 -W 3 = E'i 
E 2 

^1 
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(P, Q, i2;p, q, r, fc; M) = (1, 1, 4; 4, 2, 4, 6; 256) (3.18) 
Symmetrizing group: Sg^ = Z 2 : 

-1 1 

12r-6' y o' (3.19) 

(it, u_i,uo, ui) ^ (u- 2v , V-i - Vq, -Vq, Vi - v ) 

Polynomial invariants: 

Ei :=6u-4v_i -2u -4ui, E'i := 4ui — 2w , E 2 := 4i; 2 (3.20) 
Relation between the Halphen variables: 



(3.21) 
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6C i-> 3A: 
Rational map: 

2 8 (3G-2) 



F = G + 
f 



(G-6) 2 ' 
(9 + 3) 3 



W 2 - W 3 = 



J 3 

2E 3 



6E ^ 3B: 
Rational map: 

F = G + 



f = 



(3G + 13) 2 
(1 + 2gf 



27g< 

Powers and coefficient in (3.9): 



Symmetrization quotient: r (3)/r (6) = S\ 2 J' /S\ 2 J' = S 3 /Z 2 
Relation between the Halphen variables: 

3W! + W 2 - W 3 = Su - 5v + 3vi 
Mv - 3^) 2 



3v 1 - v 
27 



(3.22) 



270/ -1)3 
Powers and coefficient in (3.9): 

(P, Q, R- p, q, r, k; M) = (3,1, 2; 2, 3, 2, 6; 3 6 ) (3.23) 

Symmetrization quotient: Sg 2 ^ /Sg 2 ^ = S3/Z2 
Polynomial Invariants: 

Si :=6u-2u_ 3 -3i;o-2i;i, E 2 := vi(4ui - 3v ), E 3 := vl(9v - 8vi) (3.24) 

Relation between the Halphen variables: 

3Wi + W 2 + 2W 3 = Ei 

2E 2 

^ - ^ = ^ (3.25) 



(3.26) 



(P, Q, R- p, q, r, fc; M) = (3, 1, -1; 0, 5, -3, 3; — ) (3.27) 



Wl ~ W * = 9 Vl - v (3-28) 
W 2 -W 3 = ,? '" 
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9B h-> 3B: 
Rational map: 

3 3 (2G + 3) 



F = G + 
f 



G 2 + 3G + 9' 
(2 + S) 3 



9(1 + g + g 2 ) 
Powers and coefficient in (3.9): 



(P, Q, R; p, q, r, fc; M) = (3, 1, -1; 4, 4, -6, 3; 1) 

Symmetrizing group: S§| = r (3)/r (9) = Z 3 

t cog + 2 

r I— > , g I— > — , 

3r — 1 g — u 

(U, V u , Vi) l-> (li - 2l70, -Uo>, U w - Ujs, ui - Ifo) 

Polynomial invariants: 

Ei := 3-u + 6^i - 4v w - 4^, E 2 := 3(v w - i^-it, + 
E 3 :=((l-a;K-(l-a;)^) 3 

Relation between the Halphen variables: 



2Wi + W 2 = Ei 



Wi - W 3 

w 2 -w 3 



3E 3 

2a/3zE3 
3E 2 



6D ^ 2A: 
Rational map: 



_ _ 3 7 (2G 2 + 20G + 131) 
r = G H 



/ 



(G-4)3 
(4ff + 23) 4 

47(^-1)3 



Powers and coefficient in (3.9): 

(P, Q, R- p, q, r, k; M) = (2, 1, 1; 2, 2, 1, 4; 4 7 ) 

Symmetrization quotient: /SgfP = S4/S3 
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Relation between the Halphen variables: 

Wi + W 2 + W 3 = 4u - 2vf3 - 2vp - vi 

2 v / 2((-5z + v / 2)^ + (5z + v / 2)^) 

W\ — W 3 



((22a + V2)vf3 + (-22i + V2)vp) ((7i + 4>/2> j 8 + (-7i + 4y/2)vp) 
16V2vpv ((22i + y/2)v p + (-22a + V2)v $ ) 



Wo - W-x = 

((-5i + V2)v/3 + (5i + v^)^) ((7i + 4>/2> i 8 + (-7z + 4^)^) 

(3.36) 

6c 1— > 2o: 
Rational map: 

3 5 (2G 2 +81) 



F = G 
f 



G 3 

(g + 3f(g-l 



3,„ ^ ( 3 - 37 ) 



16^ 3 

Powers and coefficient in (3.9): 

(P, Q, R; p, q, r, k; M) = (1, 1, 1; 2, 0, 2, 3; 16) (3.38) 
Symmetrization quotient: Sg^' 2 /S^' 2 = A4/Z3 

(Remark: The group 18 1 2 has a function field that is not of genus 0, and hence it does not 
appear in the list of replicable functions [FMN].) 

Relation between the Halphen variables: 

W 1 + W 2 + W 3 = 3u - 2v-! - 2 Vl 

w 1 -w 3 = 



2Wo(3WO — 4:V 



\2 



w 2 -w 3 



(v -2v 1 )(3v -2v 1 ) (3.39) 
2(3wo -2 Vl ) 2 



3vq — Av\ 



6E ^ 2B: 
Rational map: 



„ _ 3 3 (10G 2 + 44G' + 43) 
F = G H 



/ = 



(G + 4) 3 
(8^2 + 2Qg _ 1)2 

(8<7+l) 3 



(3.40) 



Powers and coefficient in (3.9): 

(P, Q, R; p, q, r, k; M) = (2, 1, -1; 5, 0, -4, 2; 2 12 ) (3.41) 
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j18D' /g18D' 



'6E 1^23 



Relation between the Halphen variables: 

W 1 + W 2 + W 3 = 
W 1 -W 3 = - 

W 2 -W 3 

8E ^ 4a: 
Rational map: 



2u — hv_i + Av\ 

8 

2(v 2 + 18v vi - 27i> 2 ) 
vq — 9f i 
128ugui 

(u -9vi)(vg + 18v vi -27w 2 )' 



/ 



2 4 (5G 2 - 16) 



G(G 2 - 16) 
i(g-i) 4 



8g(g 2 - 1) 



Powers and coefficient in (3.9): 

(P, Q, R; p, q, r, k; M) = (1, 1, 2; 2, 2, 2, 4; 64) 
Symmetrizing group: Sff' = Z 4 

0-1 



-1 







8r-4' •' o + l 
1,^0,^1) !->•(«- 2v_i,v - v-\,v\ -v-i,-v-i) 

Polynomial invariants: 

Ei := Au — 2v-i — 2vq — 2vi, E2 := (i>o — V-\ — Vi) 2 
E3 := (vq - v-i - vi)(v + iv-i - ivi) 2 

Relation between the Halphen variables: 

W 1 + W 2 + W 3 = 

w 1 -w 3 = 



w 2 -w 3 = --^ 



Si 

E 2 
"E 3 

E_3 

E 2 



4. n + 1— Vertex Systems. Further Remarks. 



(3.42) 



(3.43) 



(3.44) 



(3.45) 



(3.46) 



(3.47) 



The systems listed in the preceding sections comprise a small subclass of the set of Haupt- 
moduls in [FMN]. These functions all satisfy Schwarzian equations of type (1.10) with R(f) 
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of the form (1.37), (1.38) for 1 < n < 25 and determine solutions of generalized Halphen 
equations of the type introduced in [02]. In the case of n finite poles in the rational function 
R(f), at points {cij}j = i 5 ... n , the quantities {i>i}i=i 5 ... n defined, as in eq. (3.1), by 

u:=X = l^, Vl : = ±(X - X t ) = i = l,...n, (4.1) 

satisfy the set of quadratic constraints 

(di - cij)viVj + (aj - a k )vjV k + (a k - ai)v k Vi = 0. i,j,k=l,...n, (4.2) 

which span an n — 2 dimensional space of quadratic forms vanishing on the quantities {vi}. 
They satisfy the differential equations 

v[ = -2v1 + 2uvi, i = l,...n, (4.3) 

while the remaining phase space variable u satisfies 



u' = u 2 



~ Yl r H V i V 3 ( 4 - 4 ) 



when R(f) is of the form 



1 n 

(45) 



■,3 

In (4.5) we have the freedom of adding any linear combination of the vanishing expressions 



{f -ai){f -aj) + (f -aj){f -ak) + {f ~a k )(f -di) °' ^' 6 ' 1 

which amounts to adding linear combinations of the quadratic forms of eqs. (4.2) to the 
RHS of (4.3) and (4.4). Using this freedom, we can always choose the decomposition in 
(4.5) such that the quadratic form in eq. (4.4) defined by the coefficients be tridiagonal. 
Choosing the ordering (ai, . . . , a n , oo) corresponding to a positively oriented path on the 
boudary of the fundamental region makes this unique. The differential systems are defined on 
the 3-dimensional subvariety of the space with linear coordinates (u, v±, . . . v n ) cut out by the 
quadrics in (4.2), and are determined by the set of n+1 quadratic forms appearing in eqs. (4.3), 
(4.4), defined modulo those in eq. (4.2). We note that the n quadratic forms appearing in 
eq. (4.3) are independent of the values of the parameters, while those defining the constraints 
(4.2) depend only on the locations {ai}i=i... n of the poles. The diagonal coefficients of the 
quadratic form in (4.4) are related to the angles {aj7r}j = i... n at the n finite vertices of the 
fundamental polygon by 

r u = l-al (4.7) 
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In the case of real poles the off diagonal coefficients are determined, modulo the vanishing 
quadratic forms (4.2), by the values of the accessory parameters [GS] of the corresponding 
mapping of the fundamental polygon to the upper half plane. 

As an illustration of a system associated to one of the higher level functions, consider 
the case 72e which, taken in the normalization of [FMN], may be expressed as a ratio of 
^-functions 

= 7 7 (24r)r ? (36r) 
1 77(12t)??(72t)' 1 J 

The fundamental region has 25 finite vertices, the largest number amongst the replicable 
functions. These are mapped in the /-plane to the origin / = and to the twelfth roots of 
unity times the two reciprocal radii (y/2 + 1)5 and (y/2 — 1)3. The rational function R(f) 
entering in eqs. (1.3) and (1.10) is given by 

R(f) 1 | 864/^ + 1)^-/4 + 1)2 

U ) 4/2 f (/6 + 2/3 _ l)2 (/ 6 _ 2/3 _ 1)2 (/ 12 + 6/ 6 + 1)2 ' ^ " > 

We denote the finite poles as 

(m — 1)-kx i — 1 (m—l)iri , — i 

a := 0, a m :=e a (V2 + l)s, ai 2+m := e 5 (V2-1)*, m = l,...12. 

(4.10) 

In the notation defined above, the 25 functions {vo,v m } m =i...24 satisfy the usual constraint 
equations (4.2), and equation (4.3) for the derivatives of the ^-variables. The quadratic form 
entering in the eq. (4.4) for u' is 

n g 24 g 11 

^2 r i3 V i V 3 = V + ^ V rn ~ g X] ( X ~ (( 2 + V^)u m «m+1 + (2 - V2)v 12+m V 13+m ^ . 

i,j = l m=l m=l 

(4.11) 

The two angles in the fundamental polygon at the vertices mapping to and oo therefore 
vanish, while the others are tt/2. The simplicity of this expression is due to the invariance of 
the Schwarzian derivative under the transformation r i— > r + y^, which generates the cyclic 
group action 

/^e^/, R(f)^e-^R(f), (4.12) 
and to the inversion symmetry 

R(l/f) = fR(f). (4.13) 

As a final remark we note that, for the general n finite pole case, when R(f) is of the 
form (4.5) there is an equivalent way of expressing the Schwarzian equation (1.10) and the 
associated system (4.2)-(4.4) in terms of an unconstrained dynamical system on the SL(2, C) 
group manifold. To do this, let 

9(t)'-=(c £)eSL(2,C), AD-BC=1 (4.14) 
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denote an integral curve in SL(2, C) for the equation 



where 



7 



AC 2 (Ca i + D)(Ca j + D) C 2 \ C 

Defining /(r) to be 

f •= - 



(4-17) 



it follows that this satisfies (1.10), and that 

satisfy the system (4.2)-(4.4). Equivalently, Ohyama's variables (1.39) are obtained by apply- 
ing (7 as a linear fractional transformation to {oo, a\, . . . a m } 

Conversely, up to a choice of branch in (/')^ (which does not affect the projective class of g), 
we may always express a solution / of (1.10) in this way by defining g to be 

, ' (£1 2/ ' 2 -{/"V (4.20) 



2(/')* V 2 /' -2//' 

A number of further questions suggest themselves in relation to this work. The first 
concerns the origin of the Fuchsian equation associated to each of these Hauptmoduls. For 
the case A(r) discussed in the introduction, the elliptic integral formulae (1.15a), (1.15b) show 
that the associated hypergeometric equations are Picard-Fuchs equations corresponding to the 
family of elliptic curves parametrized, in Jacobi's form, by A(r). A similar interpretation was 
found by Ohyama [03] for the 4-vertex example denoted here as 9S, for which the underlying 
Fuchsian equation is a Picard-Fuchs equation for the Hesse pencil of elliptic curves. All the 
hypergeometric functions have Euler integral representations, but it is not clear whether their 
equations might similarly be interpreted in terms of parametric families of elliptic curves. The 
same question may be asked for the other Hauptmoduls considered here. Another question 
that naturally arises is: what are the analogues of the Chazy equation that are satisfied by 
the logarithmic derivatives of the analytic forms entering in eqs. (2.7) and (3.9), and how may 
these forms be expressed in terms of some standard set such as, e.g., the modular discriminant 
appearing in eq. (1.27), Eisenstein series, or simply as ratios of ^-functions? 
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A. Appendix 

We recall here the standard definitions [A, F, WW] of the functions $2, $3, $4 and ?y, and 
a number of properties and relations between them that are helpful in verifying some of the 
formulae of Sections 2 and 3. 

Dedekind ^-function: 

00 

V(r):=q^l[(l-q n ), (q ■= (Al) 

n=l 

Null ^-functions: 

00 00 
2 (r):= <l Un+lf =2g5 JJ(l_g")(l + g») 2 (A2) 

n= — 00 n=l 

00 00 

Mr):= £ ^ n2 =II( 1 -^)( 1 + ^" i ) 2 ( A3 ) 

n= — oo n=l 

00 00 

n= — 00 n=l 

Relations between t? 2) "#3 7 ^4 7 an d r\: 

Mr) = 2^ (A5) 

V 3 (r) = \&2{r)Ur)Mr) (A8) 



Modular transformations of $2, $3, #4, and ?y: 



t? 2 (t + 1) = e^ 2 (r), ^ 2 (-l/r) = (-zr)^ 4 ( r ) (A9) 

i? 3 (r+l) = 7?4(r), i? 3 (-l/r) = (-ir)^ 3 (r) (A10) 

tf 4 ( r + 1) = tf 3 ( T ), tf 4 (-l/r) = (-it)H 2 (t) (All) 

7/(r + l) = e^?/(r), 7/(-l/r) = (-zr)^(r) (A12) 



34 



Modular solutions to generalized Halphen equations 



Identities satisfied by $2, $3, #4, V- 



^(t) = 0*(t) + 01(t) (A13) 

20§(2t) = 0§(t) + 02(t) (A14) 

#1(t) = 20 2 (2t)# 3 (2t) (A15) 

^(2r) =^ 3 (r)^ 4 (r) (A16) 

77(r)77(4r) 

^(r + 1) = e ®?7 3 (r) - 3VSe~ ^ r] 3 (9r) (A18) 



Differential relations satisfied by $2, "$3, ^4: 



tf 2 ^3 4 

_ Ha - l — 

$2 #4 4 



^ (A20) 
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